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1. Introduction 

Two-species atomic Bose-Einstein condensates (BECs) exhibit notably richer dynamical 
phenomena than single- species BECs. The inter-species interactions between the two 
components affect the nonlinear dynamics of the individual BECs and the two-species mixture 
may be in a miscible or immiscible phase [[I]|2l|3l|4l|51|6l|71[8l, exhibiting, e.g., spin [9] and 
shock waves [HOl, vector solitons [|I3[II][l2l[l3l[l4l[l5l[l6l[l7lll8land other topological 
defects and textures [[T9l . In optical lattices a bosonic two-component mixture has attracted an 
increasing experimental interest, e.g., in controlled collisions and multi-particle entanglement 
[|20]| . in mixing of ^^Rb and "^^K in 3D lattices [21fl. in super-exchange interactions [22J, in 
spin-gradient thermometry ["231, and in sub-shot-noise quantum interferometry [24 J . 

In this paper we study both numerically and analytically the effects of quantum and 
thermal fluctuations on a two-species BEC when the condensates are confined in an optical 
lattice in a highly-elongated ID trap and the lattice potential is slowly tumed up. In single- 
species bosonic atomic gases the interplay between enhanced quantum fluctuations in an 
optical lattice and the repulsive inter-atomic interactions has experimentally been shown to 
result in strongly reduced atom number fluctuations and the loss of phase coherence between 
the atoms in different lattice sites [l25l|2S|271IMll29l|30l|3l]|32l|33l. The reduction in atom 
number fluctuations has been exploited in preparation of spin-squeezed states f^Q'l that are 
suitable for quantum-enhanced interferometry. The phase separation dynamics of a two- 
species harmonically-trapped BEC was experimentally observed in two immiscible hyperfine 
components of ^^Rb [1], due to the long lifetimes of the two-fluid system that results from a 
fortuitous cancellation of the scattering lengths [jSU. Controllable spatial separation dynamics 
has recently been observed in a mixture of ^^Rb and ^^Rb atomic BECs by tuning the inter- 
species interactions with a magnetic Feshbach resonance [|2l and by introducing a linear 
electromagnetic coupling between the two internal states [|35l that creates an effective dressed 
state description for the two components [|36ll37l . 

In our simulations the atoms are initially confined in a shallow lattice and we 
continuously turn up the lattice potential. The enhanced effective interactions result from the 
reduced hopping amplitude of the atoms along the lattice and quantum fluctuations become 
more dominant in a deep lattice. We calculate the on-site atom number fluctuations in 
individual lattice sites and the inter and intra species relative phase coherence between the 
atoms in different sites. The numerical results are compared with the analytically calculated 
values that we derive in the appendix using the Bogoliubov theory. Even in the miscible 
regime of the two-species BEC system quantum fluctuations eventually destroy the long- 
range coherence of the atoms along the lattice, fragmenting the condensates. The inter and 
intra-species relative phase coherence between the atoms even in the adjacent sites is notably 
reduced in deep lattices, but the inter-species coherence remains higher close to the onset of 
the phase separation instability. We find that the repulsive inter-species interactions increase 
the inter-species relative phase coherence, but have only a weak effect on the intra-species 
coherence and the on-site atom number fluctuations. The coherence typically stabilizes to 
a non- vanishing finite value after the lattice ramping and we evaluate its spatial dependence 
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along the lattice, demonstrating a clearly reduced spatial coherence length of the system. In 
the dynamically unstable regime we find considerably enhanced atom number fluctuations, 
stronger loss of phase coherence, and the spontaneous formation of metastable configurations 
of interleaved domains of the two spatially- separated components. We calculate quantum 
mechanical expectation values and uncertainties of the domain length and find that they 
depend on the strength of quantum fluctuations, deviating from the classical mean-field values. 

The experimentally observed phase separation dynamics of Ref. [j2l in the uniform space 
was theoretically studied using the classical mean-field theory in Ref. [|38]| . A multiorbital 
wavefunction analysis was employed in Ref. ||39l to demonstrate that stronger inter-species 
interactions lead to a shorter domain length. There has been an increasing interest in the 
experimental studies of spontaneous symmetry breaking and pattern formation also in other 
ultra-cold atomic systems [i40l ITT] and the two-species condensate with a coupling between 
the two internal states [35^ has been proposed as a system to study the Kibble-Zurek defect 
formation mechanism in phase transitions B21I43J . 

In order to perform efficiently the numerical simulations in a lattice, we develop an 
approximate method to describe the non-equilibrium dynamics of a two-species condensate 
mixture that is based on the truncated Wigner approximation (TWA) [[44l l45l l46l l47l l48l 
|49l [soil of the stochastic phase space dynamics. The two-species BEC equilibrium state is 
solved within the classical Bogoliubov approximation where the excitations are evaluated 
in the tight-binding approximation of the uniform two-species lattice Hamiltonian. The 
amplitudes of the Bogoliubov phonon modes are then stochastically sampled according to 
a probability distribution given by the Wigner distributions of the ideal harmonic oscillators, 
as in the single-component TWA approaches [|47l . Each stochastic realization of the initial 
state is propagated in time according to the classical mean-field dynamics, so that individual 
stochastic trajectories represent potential outcomes of single experimental runs and quantum 
mechanical expectation values and fluctuations are calculated from the ensemble averages of 
the stochastic dynamics. One of the advantages of the approximate two-species model is its 
simplicity and the possibility to calculate analytic solutions to the initial state mode functions. 
The analytic approach to linearized excitations can also be used to calculate approximate 
ground state properties of the two-species system. We use this in Appendix to evaluate the 
intra-species relative phase coherence and the on-site atom number fluctuations. 

2. Theoretical Model 

2.1. Classical mean-field equation 

We assume that a two-species BEC is in a tightly-confined highly-elongated ID trap, so 
that any density fluctuations along the radial direction perpendicular to the trap axis can 
approximately be ignored. Along the axial direction the atoms experience an optical 
lattice potential that is deep enough so that the atoms can be described in the tight-binding 
approximation in which case one trap mode per lattice site is included in the dynamics. The 
classical mean-field model then follows our previous descriptions DH [161 and the equation 
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that governs the dynamics of two component BEC is the two-component discrete nonlinear 
Schrodinger equation (TCDNLSE), 

i^'l^ = + ^1^^) + f^xji^^!:^?^ , (1) 

where Jj (Jj > 0) and i/zn'' denote the nearest-neighbour hopping amplitudes and the 
wavefunction amplitude at the lattice site n of the atoms of species j (j = 1,2), respectively. 
The nonlinearities are given by the interaction coefficients Xjk U^k = 1,2) which are 
proportional to the onsite atom-atom interaction strengths and to the overlap of the lowest 
vibrational state wave functions (plP (the Wannier functions) of the two species in a given 
lattice site, i.e., 

Xii — d'Acp'^{f)t , (2) 

mj J 

and 



2nn ^JNiNjajk 
Xjk ^ - 



r d'n<^f{f^Mw (3) 

for j 4^ k. The inter and intra species scattering lengths are denoted by ajk and ajj, respectively, 
and Nj is the atom number of the species j. Here the reduced atomic mass p is given in terms 
of the atomic mass of the component nij as 

" = 7^- 

{nil + ^2) 

Since the number of atoms in each species is a conserved quantity, in the following we use 
the normalization 

= u (5) 

n=\ 

where L denotes the number of lattice sites. 



2.2. Parameter regimes of the simulations 

We study the quantum dynamics of a two-species BEC in an optical lattice. The initial state 
of the TWA simulations is generated by calculating the classical Bogoliubov modes whose 
amplitudes are sampled stochastically. The time-evolution for each stochastic realization then 
follows TCDNLSE of Eq. ([I]). The simulations involve a large parameter space. There are 
three nonlinearities Xn^ X12, and X\i in Eq. ([T]) and the hopping amplitudes 7^. for the two 
species may differ In addition, the two components can be moving with the different carrier 
momenta pk. The stochastic initial state fixes the atom numbers A'^i and A^2 (for given 22) 
and we can also vary the number of lattice sites L. The lattice potentials of the two components 
may be shifted with respect to each other and the two species could also experience different 
radial confinements. In the finite temperature examples we also vary the initial temperature 
T . In the following we will demonstrate how a simple analytic description for the initial state 
of the TWA simulations in terms of the Bogoliubov modes may be obtained whenever the two 
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atom currents are equal, i.e., for 71 sin(pi) = J2 sin(;?2)- In order to demonstrate some basic 
effects of the two-species quantum dynamics in a lattice, we concentrate on a simple set of 
parameter values for which the mode functions have especially compact analytic expressions. 
In all the numerical simulations we consider condensates with zero centre-of-mass momenta 
P\ = P2=Q and 

Ji = h, Xn =Xri-> Ni = N2 . (6) 

We show that the particular choice for the set of parameters is by no means necessary, 
however, and that the general formalism with the analytic initial state derivation is more 
general. For the selected parameter set we may investigate the main physical phenomena of 
the two-species lattice dynamics: the effect of the phase separation dynamics can be controlled 
by the ratio ;^'i2/;t'ii, quantum fluctuations and nonlinearity by ;^;'ii/A'^j and;(fii/7i. 

2.3. Classical Bogoliubov theory and stability analysis 



We can find steady-state solutions for the TCDNLSE ([T]) that represent propagating plane 
waves 



,l,ii) - ^ JiPin-oJit) 



where ojj is given by 



2 

^y = -2^;COs;?; + 2Ay^, 
k=i 



^jk = -r 



(7) 



(8) 



The carrier wave momenta, pj are quantized according to pj = Pj^ where Pj is an integer 
that takes value in the interval [-§, f). 

The linear stability analysis of the steady-state solution was performed in Ref. f8l using 
the classical Bogoliubov expansion. In the Bogoliubov approach the wavefunctions for each 
component in Eq. ^ are written as 



U) 



^(1 + u^J^e''^" - {vfVe-''"')e'^P'''-'^''^ 

y Li 



(9) 



In the limit of weak perturbations, the system of equations for u^^ and v'"^ may be expressed 
as an eigenvalue problem 



d 

Ij^^q = (^q^q, ^q = 



(2) 
Uq 

\ ^q J 



cr = 



(T, 




(10) 



where cr, denotes the 2x2 Pauli spin matrix. The elements of the 4x4 matrix are 
obtained from the Bogoliubov linearization procedure flSl. The quasimomenta q = ^ may 
be defined such that Q takes integer values except zero in the range [-■§, ■§). The eigenvalues 
of 9J{g correspond to the normal mode (excitation) frequencies Q.^ of the system that have 
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simple analytic expressions when two BECs have the same atomic currents ^{J\ sin(/)i) = 
J2 sin(p2))- In that case we obtain [l8l[T6l 

aq = ti+ ^ife±?3), (11) 

where 

ti = 27i sin(pi) ^m{q) , (12) 
represents a Doppler shift term of the excitation frequencies due to the superfluid flow, 

^2 = VL + V2,,, (13) 
h = yji^iq - "^Xqf + ^^^U^2,q^\2 ^2 , (14) 

are defined in terms of the single-condensate normal mode frequencies Vj^q 

= 6j^q C0Sipj)[6j^q COS(pj) + 2Ajj] , (15) 



and 



6^, = 47ysin2(^/2), (16) 



is the spectrum of an ideal, non-moving BEC. The flow is stable if the frequencies in 
Eq. ([TT]) are real for all q 0, otherwise there are small excitations in the system that grow 
exponentially in time. In these equations the normal mode frequencies of the two condensate 
species are coupled by the inter-species interactions and in the absence of the inter-species 
term, xu = 0^ we have two decoupled condensate spectra Eq. The simplest case is 

obtained when both BECs are in the normal dispersion regime with pi,p2 < n/2 and;^^^^. > 0. 
In that case the instability condition for the modes q reads flll 

Xu >XnX22 + q^LUiXu + jiXii) /2 + q^L^jJi/^- . (17) 

where jk = Jk cos pk- In the normal dispersion regime the instability first sets in for the modes 
for which \q\ is small and therefore q^L oc 1 /L. In the limit of a large lattice L ^ oo, we then 
obtain the criterion for instability 

xli^XuXii- (18) 

This stability condition is notably altered if one of the BECs exhibits anomalous dispersion 
due to superfluid flow [[HI. If the interaction strengths in Eq. ([T7]) are tuned in such a way 
that the instability is characterized by a single unstable mode, the two-component system 
can also be found in a state that is no longer dynamically stable but does not undergo a 
phase separation [16J. Instead, the two-species mixture exhibits a periodically appearing and 
disappearing vector soliton structure. 

In this work we only consider initially stationary BECs with the vanishing condensate 
momenta p\ = pi = We also assume that the hopping amplitudes for the two BECs are 
equal J\ = J2 = J and 6i_^ = 62,^ = = 47sin^(^/2). This simplifies the stability analysis. 
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The eigenvalue system for the linear stability analysis in Eq. (fT(}l ) can then be expressed as 9Ji^ 
given by 



^71,9 


-All 


Al2 


-A12 ^ 


-All 




-Al2 


A12 


Al2 


-Al2 


m,q 


-A22 


-Al2 


Al2 


-A22 





(19) 



with the definition rjjq = Ujj + Sq. The system exhibits two physical normal mode frequencies 

Q.^ ^ ^ife±?3), (20) 

For a positive definite 9J{g these are real indicating dynamical stability of the system. The 
corresponding dynamically stable eigenvectors satisfy the normalization condition ^\cr^q = 
1 . The eigenvalues -Q.q of Eq. (flOl) represent unphysical solutions with the corresponding 
eigenvectors satisfying the negative normalization ^Icr^q = -1. 

The BEC system becomes dynamically unstable when the normal mode frequencies 
exhibit nonvanishing imaginary parts, indicating perturbations that grow exponentially in 
time. The corresponding eigenvectors satisfy ^lo-^q = 0. The rate at which the instability 
sets in depends on the magnitude of the imaginary part of the eigenfrequency. 

We can solve the eigenvectors of Eq. (fTOl ) analytically. The expressions for the mode 
functions notably simplify when we consider the case;^fii = X22- We then obtain (for;^i2 0) 



^ = -^e,(6, + 2Aii)±26,Ai2 (21) 

,(1) ^ 4Ai2Q- + (^-)^-(Q-)^ ^^^^ 
4^2Ai2Q|[(Q;)2-(ri-)2] 
4A17Q* + (Qrf - (Q+)2 

(1) _ ^'^I2^''q ^ V^'^g^ y^^q-l 

^q,± — I K-^^J 



4^2Ai2Q|[(Q;)2-(^^-)2] 



and ufl = +u'ql, vfl = +v%. Here we have assumed, for the notational simplicity, that 
A12 > 0. In TWA simulation we generate the initial state noise for the configuration that is 
dynamically stable (we specifically consider thermal equilibrium states). For such states the 
normal mode frequencies are all real and the eigenmodes satisfy the normalization condition 

^Icr^q = [u%f - [v%f + [uf^' - [v^]^ = 1 . (24) 

The two-species BEC normal modes describe the dynamics of mixing between the two 
components as well as excitations of the total density in the system. In nonlinear regime 
we have Xjj ^ J- If we also have ~ XnXii, one of the frequencies approaches 
zero corresponding to the phase separation instability Eq. (fTSl) . We obtain in that case 
^q+ - ^Iq + ~ ^\ q-'^\ q' "^^^ cncrgy cxcltatlons then correspond to the 

mixing of the two species with only a weak variation in the total density of the two-species 
condensate. 
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In both stable and unstable regimes of the two- species mixture we can investigate 
the degree of overlap between the two species. We define the overlap integral of the 
wavefunctions for the mixture as 

2 



K{t) = 



(25) 



In the stable regime Eq. (|25l) describes the spin excitations of the two-component system. 
When the two-species interaction strengths satisfy the condition Eq. (fTVl ). the two-species 
system is dynamically unstable and undergoes phase separation, resulting in strongly reduced 
overlap integral values. A measure that can particularly well identify the phase-separation of 
the densities of the two species may be calculated from the sum 

T{t) = Y,\^'n\t)\^\^?(t)\\ (26) 



2.4. Truncated Wigner approximation 

In the TWA simulations we calculate ensemble averages of stochastic trajectories for which 
the time evolution follows the classical mean-field theory, but in each realization the initial 
state is stochastically sampled from a Wigner distribution that approximately synthesizes the 
quantum statistical correlations of the initial state. Approaches introduced in the TWA initial 
state generation in single-component BECs involve evaluating the initial state correlations 
within the Bogoliubov approximation ["ST, '471 or by solving the ground state and the excited 
state populations self-consistently within the Hartree-Fock-Bogoliubov approximation 

In order to implement the TWA phase-space model in a two-component BEC system 
we similarly assume that the two component stochastic fields 0'*^^-' obeys the classical field 
equations similar to the TCDNLSE (Eq. [U) 

i^/i^ = -Ji^fi, + ^1^,) + j^Xj0!:Wi' . (27) 

For the stochastic initial state generation we introduce an approximate model based on 
the classical Bogoliubov theory, described in the previous section, that provides simple 
analytically solvable mode functions. We write the both components j = 1, 2 as 

VAO^i^-) = cl^J'c^i^ + 6ll/^ , (28) 

where (p^^ denotes the normalized ground state solution of the BEC component j and the 
excited- state fluctuations are given by 

Here the eigenmodes {q, +) correspond to the eigenfrequencies of Eq. The mode 

amplitudes CQ ^ag^), are stochastically sampled from the Wigner distribution of harmonic 
oscillators as explained below. We consider a two-species system with equal populations 
Ni = N2 = N, with the interaction strengths satisfying = xii- In that case the modes Ug^ 
and Vqj, are given by Eqs. ([22]) and (|23l) . 
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Using the field decomposition Eq. and the mode functions, Eqs. (|^ and (1^ . 
we then stochastically sample the amplitudes of the mode functions by treating them as 
ideal harmonic oscillators whose distributions are determined by the corresponding Gaussian 
Wigner function \52i 

W{a%[a%r) = ^tanh(^,^)exp[-2|a^iptanh(4*)] , (30) 

where = Qfj/lkeT. The stochastic mode function amplitudes a''Jl produce the ensemble 
averages 

([aZT^%>w = nq,± + ^, (31) 



«'/,± = — 7 ' (32) 



where 

1 

oxpiQpicBT) - 1 ' 

is the usual Bose-Einstein distribution function. The factor 1/2 in Eq. (|3TI) results from 
the Wigner distribution that returns symmetrically ordered expectation values, providing the 
vacuum noise in each mode. For each stochastic realization the total number of excited-state 
atoms varies according to 



(33) 



q^O,v=: 

with the average number given by 

(Ni^^)= + l^i!'vl')%v + Iv^II'] (34) 

qi=Q,v=± 

The ground state amplitudes otq^ fluctuate in each stochastic realization [|50l l53l. The ground- 
state atom number is then obtained from the fixed total atom number A'^ in each atomic species, 

so that in each stochastic realization A^^'''' = N - N^-'^ and we set a^J-' = sJNc^^ + 1 /2. 

3. Numerical results 

3.1. Turning up the optical lattice 

We solve the stochastic dynamics during the turning up of the lattice potential for given initial 
conditions. We study the response of the system to the ramping so that both the tunneling 
coefficients and nonlinearity are time-dependent. We assume that at all times Ji = Jj = J and 
Xn = X22- For simplicity, we consider a situation where the both species have an equal mass 
m, so that the recoil frequencies cor are equal and are given by 

2ma^ 

where d denotes the lattice spacing. In a deep lattice the hopping amplitude 7 is then 
approximately given by [|54l 

4g-2VJ^3/4 

J = — ojR . (36) 
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Here s denotes the lattice height in the units of the lattice photon recoil energy. In a tightly- 
confined elongated ID trap the atoms are assumed to be confined in the radial vibrational 
ground state. If the radial confinement is the same for the both species, we obtain 



where denotes the trapping frequency of the radial confinement and Is = {%lm0.s)^l^, 
where Q.s - Is^^^ojr is the axial trap frequency at the lattice site minimum. 

In our simulations we numerically solve the time-evolution using the split-step 
method . As the lattice is turned up the hopping amplitude J rapidly decreases according 
to Eq. (l36l) and the interaction strength ;^f,y slowly increases according to Eq. (|37| ). due to the 
stronger confinement of atoms in individual sites. We increase the lattice height linearly at 
the rate 6, so that the lattice height satisfies s(t) = 5, -I- 6t, where 5, denotes the initial height. 
We choose 6 = 2x IO'^cor and 5, = 2, resulting in the initial value of 7 ^ O.Hojr. Unless it 
is stated explicitly, the length of the lattice is L = 64 and the number of atoms in each species 
in each run is taken to be N/L = 40. In most cases we choose the initial value for the intra- 
species interaction strengths xn = X22 - 0.60ojr, resulting in Q.j_ajj/d - 7.9 x IO'^ojr. 
The initial state fluctuations are evaluated within the Bogoliubov approximation and the 
interaction strengths are selected in such a way that the excited state population remains low. 
For;^fn = X22 - 0.60cc)r and^i2/;^fii = 0.1 at T = 0, the number of atoms initially depleted 
from the ground state due to quantum fluctuations {N^}^) = (M^^) - 70 (corresponding to 
2.7% depleted fraction). All the presented simulation results are for T = 0, except the finite 
temperature cases studied in Sec. 13.2.41 

We study the effect of ramping of the lattice both in the stable (xi2 < yjXnXii) and 
unstable ixn <: y/XnXiD regime and write the inter-species interaction strength (for i 4^ j) 
as 



The nonlinearity corresponding to the inter-species interactions is tuned by varying the 
parameter y in Eq. (l38l) . Here y < 1 corresponds to the stable regime while 7 > 1 implies the 
dynamical phase-separation instability. In cases where we study the two-species system in the 
unstable phase separation regime, we use the initial two-species mixture that is dynamically 
stable, but change the value of y from the stable to unstable regime immediately after the 
lattice ramp. 

For a fixed nonlinearity quantum fluctuations are enhanced by reducing the atom number 
(and correspondingly increasing the scattering length). For large atom numbers and in shallow 
optical lattices quantum effects are weak and the system can be accurately described by 
classical mean-field theory. As the lattice is turned up, the effective interactions become 
stronger and quantum fluctuations are enhanced. 

The ramping of the lattice is adiabatic if the rate of change of the parameters in the 
Hamiltonian is slow compared to the lowest collective excitation frequency [|56l |47]| . The 
fastest varying parameter during the ramping is the hopping amplitude and we require for an 




(37) 



Xij = y^lXiiXjj- 



(38) 
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adiabatic ramp 



1 djjt) 

Jit) dt 



< min[Q^(0] . 



(39) 



If this condition is not satisfied, the system can be excited from its ground state during the 
tuming up process of the lattice. 

3.1.1. Validity of the parameter regimes Our chosen set of parameters, explained in Sec. 12.21 
captures the essential features of the condensate fragmentation, reduced atom number 
fluctuations and the domain formation. For the phase separation dynamics, the important 
condition is that of the dynamical instability ;\fj2 > X 11X21 and the precise r:dXio xnlxii is less 
relevant. Experimentally, the interaction strengths xtj can be controlled using two-species 
Feshbach resonances [El or by introducing a linear electromagnetic coupling between the two 
intemal states [|35l . The interaction parameter ;(f,y incorporates the atom numbers and the ratio 
XnlXii may also be tuned by changing the relative atom population of the two condensate 
components. The intra-species interaction strength 12 can be controlled in a spin-dependent 
optical lattice fW\ by changing the relative lattice positions of the two species and therefore 
modifying the spatial overlap integral between the lattice site wavef unctions. 

The two-species condensate system may be realized by using two different hyperfine 
levels of the same atom 13 or by trapping two entirely different atoms, e.g., a "^^K-^^Rb 
mixture [21]. In the case of far-detuned optical lattice, the potential experienced by the 
atoms in two different hyperfine levels of the same atom is typically the same, resulting in 
the identical values for the hopping amplitudes. 

The atom dynamics can be described by a ID model if the frequency of the radial 
trapping potential Q.^ is larger than the chemical potential of the atoms ojj in Eq. ([8]) and 
the thermal energy IcbT . The typical values used in the numerics xn = X21 - 0.60a»« at 
s = 2, corresponding io xn - at 5 = 40 and Q.j_an/d - 7.9 x \0~^a>R. We therefore 

obtain for the requirement of the ID dynamic description Q.j^ > I.Sojr/L and d/an > 250. 
The first expression yields in terms of the lattice spacing d > 0.3l± where l± = ^lfil{mQ.^^ 
denotes the radial width of the lattice site mode function. For the scattering length an ^ 5nm 
we obtain from the second inequality d > 1.3/im. In ID lattices the lattice spacing can be 
easily controlled by adjusting the angle between the counter-propagating lasers that form the 
standing- wave pattem, so that d = n/[k sin(0/2)] for the intersection angle 6 and wavenumber 
k. In the simulations we assume 40 atoms per site and we require that the atom density is 
sufficiently low so that the inelastic atom losses remain weak. The three-body loss rate of the 
atoms at the site / may be approximated by 



where ^3 denotes the three-body recombination rate. In order to have a weak three-body 
loss rate n^Y/ojR «; 1, the radial trap frequency should be sufficiently weak. The parameters 
depend on the particular atom and the hyperfine state. For instance, for ^^Rb in the \F = 
l,mf = -1) state we have Kj, ^ 5.8 x 10"^° cmVs ll57l and we obtain for the condition of the 




(40) 
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three-body loss rate to be weak Q.^ «; 1.2 x 10^/s. This condition can be satisfied when the 
system dynamics is strictly ID. It should be also noted, however, that even in elongated traps 
that are not tightly confined, ID numerical model can provide a good qualitative description 
of the atom dynamics in the lattice [33 J . 

In the tight-binding approximation we assume that only the lowest energy band is 
occupied and one mode function per lattice site is sufficient to represent the dynamics. This 
approximation becomes better in deeper lattices and provides a reasonable description for 
s > 2 [|54l . We also require that the nonlinearity is smaller than the energy gap between the 
lowest two energy bands ~ %Q.s = 2^/s'hoJR, so that the higher bands are not occupied. This 
yields 2 yfscoR > Xtj/L which is well satisfied for all the studied lattice heights. 

Experimentally, the atoms are trapped in a combined optical lattice and a harmonic trap. 
The harmonic potential introduces a non-uniform atom density. This influences the phase 
separation dynamics; the condensate component with a weaker nonlinearity energetically 
favours higher density regions close to the centre of the trap [fTl|. The phase coherence and 
number fluctuations in a harmonic trap depend on the spatial location with quantum and 
thermal fluctuations stronger close to the edge of the atom cloud [l33l . Some other possible 
effects on the phase coherence are addressed in Sec. 13.2.11 A lattice with a uniform density 
and periodic boundary conditions may be realized in a toroidal trap with an optical lattice 
formed by the interference of two counter-rotating Laguerre-Gaussian laser beams [|58l . 

3.2. Dynamically stable regime 

We first consider the two-species BEC dynamics in an optical lattice in the dynamically stable 
regime of the spatially overlapping condensate mixture for y < 1. This corresponds to the 
situation where the inter-species interaction is not strong enough to cause the phase separation 
of the two components and all the normal mode frequencies in Eq. (fTTI) are real. 

3.2.1. Condensate fragmentation and phase coherence The atoms are initially confined in 
a shallow lattice and we continuously turn up the lattice potential. The effect of quantum 
fluctuations on atom dynamics in the lattice can be studied by calculating the phase coherence 
between the atoms in different lattice sites. When the lattice is turned up the hopping 
amplitude of atoms between adjacent sites rapidly decreases, resulting in the reduction of 
kinetic energy of the atoms and hence stronger effective interactions. Quantum fluctuations 
in the system are enhanced and the phase coherence between the atoms in different sites is 
destroyed as the condensates undergo fragmentation. We evaluate the loss of phase coherence 
by calculating the absolute value of the normalized relative intra-species phase coherence 
between the atoms in different sites, 

C,^i = -P=^, (7 = 1,2), (41) 

between the atoms of the same species j in sites k and /. Since we choose N\ = N2 = N and 
Xn = X22, the two species have identical nonlinear properties that are spontaneously broken 
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Figure 1. Relative phase coherence between the atoms in the same atomic species at two 
adjacent sites C\ (left) during and after the turning up of the lattice potential for different final 
lattice heights i = 10, 20, 30, 40 (curves from top to bottom) with inter-species nonlinear 
parameter y = 0.1 (top row) and y - 0.5 (bottom row). In all the cases the lattice is turned 
up at the same rate and the end of the ramping times are marked on the x axis using the 
same colours as in the corresponding curves. The corresponding stationary (averaged) values 
of the coherence C\ (right) as a function of lattice height. The coherence C\ monotonically 
decreases as the lattice is ramped up. The initial lattice height is i, = 2 the corresponding 
nonUnear interaction parameter is 11 = 0.60wr. 



only due to nonlinear interactions, e.g., in the phase separation. We describe the relative phase 
coherence between the atoms in the two different species by 

^a2) ^ (42) 

In Fig. [H we show the relative intra-species phase coherence between the atoms in two 
adjacent sites C\ during and after the turning up of the lattice potential. The different curves 
correspond to different values of the final lattice height 5= 10, 20, 30 and 40. The displayed 
cases have the inter-species interaction strength y = 0.1 and y = 0.5, defined by Eq. (1381) . 
We also show the corresponding stationary (averaged) values of the coherence Ci that are 
obtained after the turning up of the lattice potential. These demonstrate how the intra-species 
coherence rapidly decreases as the lattice becomes deeper, indicating an increasing degree of 
fragmentation of the BEC as a function of the final lattice height. We only find a very weak 
dependence of the intra-species relative phase coherence Ci on the inter-species interaction 
strength y for the values of y in the stable regime we considered in the simulations (from 
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Figure 2. Relative phase coherence between the atoms separated by 10 lattice sites and 
belonging to the same atomic species C\q (left) for different final lattice heights s = 
10, 20, 30, 40 (curves from top to bottom) with inter-species nonlinear parameter y - 0.5. 
The end of the ramping times are marked on the x axis. The corresponding stationary 
(averaged) values of Cio (middle) as a function of lattice height. The spatial dependence of the 
relative phase coherence along the lattice between the atoms in the same species, displaying 
the stationary (averaged) values of Cj as a function of the relative lattice position j (right). We 
specifically show Cj for j - 1 , 5, 10, 15, 20, 25, 30 for different cases of the final lattice height 
s - 10,20,30,40 (curves from top to bottom). The coherence Cm decreases notably more 
rapidly than the coherence C\ between the atoms in adjacent sites as a function of the lattice 
height. As in the case of C\, we find that C\q depends only weakly on y in the stable regime 
for the values we studied from y = 0.1 to 0.75. 



7 = 0.1 to 0.95). In Ref. [ j2T1 | the presence of "^'K was found to lead to lower visibility of 
the interference fringes of the ^''Rb in the two-species mixture in an optical lattice. In the 
experiment, however, the interactions increased the atom density of ^^Rb close to the trap 
centre due to the inhomogeneous trapping potential and drove the system closer to the onset 
of the Mott insulator transition, as also demonstrated in Ref. Il59i where the adiabaticity and 
the visibility of interference fringes in loading of a bosonic mixture to an optical lattice was 
studied using Gutzwiller mean-field method. 

The intra- species long-range spatial coherence is shown in Fig.|2l We display the relative 
phase coherence Cio between the atoms in one of the sites and in its 10th nearest neighbour 
site. This decays notably faster than the coherence C\ between the atoms in the adjacent sites 
in Fig. [U The spatial dependence of the relative phase coherence along the lattice is also 
shown in Fig. [2l We calculated the stationary, averaged values of Cj for different j when the 
value of the coherence was stabilized after the end of the lattice ramp. The graphs show the 
decay of the spatial coherence along the lattice. The coherence exhibits a very slow decay for 
large values of j and remains high for the case of shallow lattices. 

The numerically calculated values of the intra-species relative phase coherence may 
be compared to the analytic estimates obtained for the ground state of the optical lattice 
system in Appendix. The results for the nearest-neighbour phase coherence for the linearized 
Bogoliubov theory of the fluctuations in the ground state, displayed in Fig. lAll in Appendix, 
provide a good agreement with those obtained in the numerical TWA simulations, shown in 
Fig. [U The long-range coherence along the lattice in the TWA numerics, however, decays 
more rapidly as a function of the spatial separation than in the case of the ground- state 
calculation. 
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Figure 3. Relative phase coherence between the atoms in the different atomic species at two 
adjacent sites C*'^' for different final lattice heights s = 10, 20, 30, 40 (curves from top to 
bottom) in the stable regime with inter-species nonlinear parameter y = 0.1, 0.25, 0.5 and 0.75. 
The other parameters are the same as in Fig.[T] 



The inter-species coherence C\ is shown in Fig. [3]for different final lattice heights and 
for different values of inter-species interactions, y = 0.1, 0.25, 0.5 and 0.75. Although the 
intra-species coherence Ci is not strongly affected by the inter-species interaction strength y, 
the inter-species coherence C\^^'' is very sensitive to y even when the two-species mixture is 
miscible. In particular, the relative inter-species coherence Cj'^^ is enhanced due to the inter- 
species interactions whenxn is increased. C[^^^ becomes high as the system approaches to the 
onset of the phase separation instability, as shown in stationary averaged values of Fig. |4]that 
are calculated after the turning up of the lattice. This is because the effective interactions in 
a perfectly overlapping two-species mixture are almost completely canceled out immediately 
below the onset of the instability for;t'n ^ X22 - Xn- We also show the decay of the spatial 
coherence along the lattice by displaying C^.^^^ for different values of j in Fig.|4](on right). The 
stationary, averaged values for the coherence are calculated after the end of the ramp when 
the coherence has stabilized. 

3.2.2. Atom number fluctuations The stochastic simulations in the Wigner representation 
yield expectation values for the operators that are symmetrically ordered. These can be 
transformed to normally ordered expectation values of the atoms in each lattice site. For 
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Figure 4. The stationary averaged values of the relative inter-species phase coherence Cj'^^ as 
a function of the inter-species interaction strength y (left) for different values of the final lattice 
height s = 10, 20, 30, 40 (curves from top to bottom) and as a function of the lattice height 
(middle) for different values of the inter-species interaction strength y - 0.1,0.25,0.5,0.75. 
The stationary values are obtained after the phase coherence is stabilized after the end of the 
ramping. C*'^^ increases as the inter-species interaction strength increases in the miscible 
regime. The spatial dependence of the inter-species relative phase coherence along the lattice, 

(12) 

displaying the stationary (averaged) values of Cj as a function of the relative lattice position 

(12) 

j (right). We show C ^ for j - 1,5, 10, 15,20 for different cases of the final lattice height 
s = 10, 20, 30, 40 (curves from top to bottom). 



instance, the atom number expectation values and fluctuations are obtained from 

{nf) = {[iljfr^f)-\, (43) 



for each species- j at the lattice site I. 

It is useful to scale the atom number fluctuations to those obtained in the Poissonian limit 
that correspond to the fluctuations resulting in an instantaneous splitting or in the splitting of 
a non-interacting gas 

An^j^ 

'sqz 



nl^U^^, (45) 



so that the values n^qz < 1 indicate reduced on-site atom number fluctuations. The on- 
site and the relative atom number fluctuations between the atoms in different lattice sites 
were calculated within TWA for a single- species BEC and compared with experimental 
observations in Ref. [33] providing a good qualitative agreement. 

In Fig. [5] we show the scaled on-site atom number fluctuations ^sq^ in one of the sites for 
the different interaction strengths. The stationary averaged values of atom number fluctuations 
are shown in Fig. |6l We show both the dependence of the fluctuations on the lattice height as 
well as on the inter-species interactions. The atom number fluctuations are strongly reduced 
as the final lattice height is increased. Strong suppression on number fluctuations for deep 
lattices is associated with the enhanced phase fluctuations and the condensate fragmentation, 
and the reduced atom number fluctuations correlate with the previously calculated values of 
the loss of phase coherence. We also find that the inter-species interactions generally enhance 
the atom number fluctuations. 



A two-species Bose-Einstein condensate in an optical lattice 



17 



N 
IT 
V> 




N 

m 



40000 




40000 



Figure 5. Scaled on-site atom number fluctuations Msqz in one of the lattice sites in the stable 
regime for diff'erent final lattice heights s - 10,20,30,40 during and after the turning up of 
the lattice for different values of the inter-species nonlinearity y = 0.1 and 0.75. The other 
parameters are same as those in Fig. [T] Due to the symmetry x\\ - X22 the average atom 
number fluctuations are the same for the both species of the two-component system. 
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Figure 6. The stationary averaged values of the scaled on-site atom number fluctuations 



Msqz as a function of lattice heights s (left) for dififerent values of inter-species nonlinearity 
y - 0.1,0.75 and as a function of the inter-species interaction strength 7 (right) for different 
lattice final heights s - 10, 20, 30, 40 (curves from top to bottom). The parameters are the 
same as in Fig.|5] The number fluctuations are strongly reduced as a function of the final lattice 
height, but the dependence of the number fluctuations on the inter-species interactions is less 
significant. 



The non-equilibrium dynamics of the TWA simulations may again be compared to the 
analytic estimates obtained for the ground state of the optical lattice system in Appendix. The 
agreement between the linearized ground- state results of Fig. lAll in Appendix and the TWA 
results of Figs.|5]and|6]is very good, indicating that the effects of nonadiabaticity in the lattice 
ramping on the on-site atom number fluctuations are weak. 

3.2.3. Adiabaticity and excitations of relative atom populations In Fig. |7] we show the 
population of the first five lowest momentum Fourier modes for the inter- species interaction 
strength 7 = 0.1 as a function of time and a snapshot momentum distribution of all the modes. 
The four figures correspond to the value of lattice height = 10, 20, 30, 40. Turning up of 
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Figure 7. The population of first five lowest momentum Fourier modes during the turning up 
of the lattice potential for diff'erent final lattice heights s = 10, 20, 30, 40 (top two rows). Here 
y - 0.1. The lowest mode is initially occupied while the population of the higher modes is 
negligible. In the case of a deep lattice the population of excited modes is increased, indicating 
nonadiabatic turning up of the lattice potential. The ensemble-averaged populations of all the 
momentum modes at f = 20000/(yR (bottom row). Only the lowest momentum modes are 
occupied. 



the lattice potential results in a decrease in the lowest mode population. The effect is stronger 
when the lattice becomes deeper, indicating breakdown of the adiabaticity in turning up of 
the lattice. Consequently, the deeper the lattice, the larger the depletion in the lowest mode. 
The breakdown of the adiabaticity in deep lattices may be understood from Eq. (|39] ) since the 
frequency of the lowest phonon mode is reduced as the lattice becomes deeper and it becomes 
progressively more difficult to turn up the lattice adiabatically. 

We may estimate the adiabaticity of the turning up of the lattice potential using the 
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expression Eq. In the TWA simulations the lattice is turned up at the rate s{t) = si + 6t, 
with Si = 2. For y = 0.1 and 5 = 2.0 x \0~^a)R we obtain D}f\t)l^{t) ^ 10 at 5 = 2 and 0.7 
at the end of the deep lattice ramp s = 40. The condition - 1 is reached at about 

s - 34. The adiabaticity condition is much easier to violate close to the onset of the phase 
separation. For y = 0.95 and 6 = 2.0 x IO'^cor we obtain - 4- at s = 2, reducing 

to 1 at 5 - 12, and to 0.2 at ^ = 40. 

In Fig. [8] we show the the relative phase coherence between the atoms in the adjacent 
lattice sites and the on-site atom number fluctuations for three different speeds of the lattice 
ramp, representing 6 = 2.0xl0~^a»s, l.OxlO^^a*^ and 0.67 x 10~^a»«. Fory = 0.1 and the two 
slowest ramp cases the condition Q™°(?)/^(0 - 1 is never reached during the turning up of 
the lattice potential. For 0.67 x \0~^ojr the initial value - 29 at 5 = 2 is reduced 

to about 2.1 at 5 = 40 for y = 0.1 and from about 12 at 5 = 2 to about 0.5 at 5 = 40 for 
y = 0.95 - 1 is reached at about s - 29). 

Despite the improvement in the adiabaticity condition between the three difl'erent ramps, 
there are very little changes in Fig. |8] especially in atom number fluctuations. Reaching the 
limit where the ratio Q™'°(0/^(0 » 1 is both numerically and experimentally demanding 
in the case of optical lattices with large occupation numbers and a large number of sites. 
For instance, with the present parameter values, maintaining Q.^^"(t)/^(t) > 10 during the 
entire ramp to * = 40 for y = 0.95 would already require a very slow ramp speed of 
6 - 0.2 X lO'^ojR. The difficulty of adiabatically tuming up a lattice potential in the 
experiments can severely limit possibilities to reach the superfluid Mott-insulator transition 
in the case of large occupation numbers [|47l . Achieving a strong reduction in atom number 
fluctuations in the case of many atoms has been experimentally challenging even in the lattice 
systems of only a few sites [|29l|. 

The excitation of higher modes [Fig. |7l induced during the ramping process indicate 
a nonadiabatic tuming up of the lattice. The lattice-induced excitations are also reflected in 
relative atom populations between the two species (spin- 1/2 waves). We show such excitations 
in the stable regime y < 1 in Fig. [9] by displaying the overlap integral, defined by Eq. (|25l) . 
between the wave functions of two atomic species for various values of inter-species nonlinear 
parameter y = 0.1, 0.25, 0.5, 0.75. In deeper lattices the population difference is clearly 
increased. The effect of interactions, however, is again reduced as y ^ 1. 

3.2.4. Effects of temperature Finite temperature in non-equilibrium quantum dynamics 
introduces additional noise in the system, increasing the atom number fluctuations. In 
experiments on atom number squeezing and reduced on-site atom number fluctuations in 
optical lattices the finite temperature has been an important factor limiting the achievable spin 
squeezing and the suppression of the atom number fluctuations ||29ll33]| . Here we demonstrate 
the effects of the initial temperature of the atoms on the coherence properties of the two- 
species system as the lattice potential is tumed up. The temperature can be incorporated in the 
stochastic sampling of mode populations according to Eq. (|32l) in which case the width of the 
Gaussian stochastic distribution for the sampling of the initial state is increased due to thermal 
population of each phonon mode. In Fig. [TO] we show the effect of variation of temperature 
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Figure 8. The relative phase coherence between the atoms in the adjacent lattice sites Ci and 
the on-site atom number fluctuations for three different lattice ramp speeds 5 - 2.0 x IO^^ujr, 
1.0 X IO^^ur and 0.67 x lO^^ws. The lattice is turned up according to sit) - Si + 6t. We show 
two diff'erent cases of y = 0.1 (top row) and 0.95 (bottom row). 

on the relative intra-species phase coherence between the atoms in different lattice sites and 
on the on-site atom number fluctuations. The inter-species interaction parameter y = 0.5. The 
corresponding stationary averaged values of atom number fluctuations ngq^ and the relative 
intra-species phase coherence between the atoms in the adjacent lattice sites Ci are shown in 

Fig.im 

3.3. Unstable regime 

In the previous section we considered a two-species condensate mixture in the regime where 
the spatial overlap of the two species is dynamically stable, corresponding to the values 
of the inter-species interaction strength y < 1. When the parameter y is increased the 
system becomes dynamically unstable and the normal mode frequencies of Eq. (fTT)) exhibit 
nonvanishing imaginary parts, indicating perturbations that grow exponentially in time. The 
instability criteria in different regimes for static and moving condensates were analyzed in 
detail in Refs. |f8l fT6| and also the effects of matter- wave grating of the other species [[60l 
have been studied. Phase separation is a generic phenomenon that occurs in different forms 
of matter. The phase separation instability condition for a two-species BEC in a lattice is 
analogous to the phase separation instability criterion of the two BEC components that occurs 
in free space when the square of the inter-species interaction coefficient exceeds the product 
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Figure 9. Overlap integral between the two condensate species [as defined by Eq. ( |25]) 1 in the 
stable regime for different values of inter-species interactions y - 0.1,0.25,0.5,0.75 for the 
final values of the lattice height s - 10, 20, 30, 40 (curves from top to bottom). 




Figure 10. Effect of temperature on the on-site atom number fluctuations ngq^ and on the 
relative intra-species phase coherence between the atoms in the adjacent lattice sites C\. 
Here the nonlinearities xn - X22 - 0.6(yR, the lattice height s=30 and y - 0.50. The 
number of excited-state atoms per condensate component for ksT ftuDR - 0.01 and 0.02 are 
approximately 370 and 870, respectively. 
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Figure 11. The stationary averaged values of the relative intra-species phase coherence 
between the atoms in the adjacent lattice sites C\ and the on-site atom number fluctuations 
Hsqz a function of temperature, corresponding to Fig. [10] 

of the intra-species interaction coefficients. 

Here we consider the unstable regime of (y > 1) by first evaluating the thermal 
equilibrium state of the atoms in the initial state for some value of y < 1 corresponding to 
a stable regime of overlapping two-species mixture. Stable initial ground-state configuration 
allows us to evaluate the statistical noise for the initial state of the TWA simulations within the 
Bogoliubov approximation. We then continuously turn up the optical lattice potential as in the 
dynamically stable case and immediately after the final lattice height is reached we change the 
inter-species interaction y to the unstable regime. Varying the final lattice height and the atom 
number provides then information about the dependence of the phase separation dynamics on 
the lattice parameters and on quantum fluctuations. Experimentally the manipulation of the 
scattering lengths for BECs in order to drive the system from the stable to unstable regime has 
been realized for a ^^Rb-^^Rb condensate mixture using a Feshbach resonance [2J. Related 
experiments on single-component BECs by rapidly changing the scattering length from the 
stable positive to unstable negative value have generated a condensate collapse and the 
formation of bright solitons [I6ni62ll63l . In two-species condensates the efi'ective interaction 
strengths have been manipulated between stable and unstable regimes using dressed atomic 
states by electromagnetically-induced Raman transitions between the internal atomic states 

3.3.1. Domain formation The system develops instability when the inter-species interaction 
exceeds a threshold y ^ 1 determined by the intra-species interactions. After the interaction 
parameters are switched to the unstable regime at the end of the lattice ramping the atom 
densities of the two BECs show violent phase separation dynamics and individual sites 
become dominantly occupied by single species alone. In the ground-state configurations of 
phase-separated systems one species typically forms a shell around the other one, minimizing 
the surface area between the two components [|38l . In an optical lattice system we consider, 
quantum noise and the nonlinear interactions spontaneously break the symmetry of the 
uniform spatial configurations of the two species. Individual stochastic realizations of the 
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Figure 12. Representative individual stochastic realizations of the atom density distributions 
for the two-species condensates in the dynamically unstable phase separation regime that 
represent possible outcomes of single experimental runs. The two curves correspond to the 
densities of the two atomic species. The different figures correspond to different final lattice 
heights s - 5, 10, 15, 20 after the system has reached a metastable density configuration. 
The initial value of the interspecies interaction parameter during the turning up of the lattice 
potential is in the stable regime with y - 0.95. After the ramping up of the potential it is 
changed to the unstable value y - 1.1. The unstable dynamics results in the spatial phase 
separation pattern of interleaved density domains. Here the lattice size L - 128. 



TWA dynamics that represent possible outcomes of individual experimental runs show density 
domain formations of the two species. The system settles down to a metastable configuration 
of several interleaved spatially- separated density domains of the two species where the entire 
system may consist of multiple domain boundaries along the lattice. Individual sites are 
typically dominantly occupied by one atomic species alone, except close to the domain 
boundaries for the case of large domain length. In Fig. [T2]we show typical single realization 
results for the atom density distribution for lattice heights s = 5, 10, 15, 20 after the system 
has phase-separated. The number of domains increases as the lattice height increases. 

The observed phenomenon is related to the experiments on a harmonically trapped ^^Rb- 
^^Rb condensate mixture in the absence of a lattice [2] in which case one of the species 
was found to split into multiple separated atom cloudlets that appeared as distinct holes in 
the density distribution of the other species. The condensation experiment was theoretically 
analyzed in Ref. |[38l and it was demonstrated how the dynamics leads to continuous 
separation of the two species into smaller and smaller domains. As shown in Ref. [fT6ll . 
the formation of only a few density spikes and holes in the phase separation dynamics can 
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be identified as a spontaneous generation of bound pairs of dark and bright solitons Wl\ . 
The phase separation provides a mechanism for the background of a dark soliton (density 
hole) in one species to stabilize a bright soliton (density spike) in the other species due to 
an effective trap that results from the repulsion between the two species [13J. In the present 
system the stabilization of the small domains is similar to the stabilization in energetically 
metastable particle-like solitons [13] and results in metastable configurations that are not 
necessarily energetically close to the ground state. The formation of the metastable states is 
a non-equilibrium process and the interleaved pattern does not represent a thermal state. One 
should compare the observed phase- separated state to the ground state that is a maximally 
phase- separated state and minimizes the surface area between the two components by forcing 
them to the opposite sides of the trap (one component to the right and the other one to the 
left). 

In the spontaneous pattern formation, due to an instantaneous switch of the interaction 
strengths, the domain length is expected to be approximately determined by the wavelength of 
the phonon mode with the largest imaginary part of the eigenfrequency, since this corresponds 
to the unstable eigenmode that grows most rapidly (provided that the perturbation is strong 
enough to populate this mode). We can calculate the wavenumber of the mode ^max 
analytically from the expressions of Q.q. The value of ^^ax represents the fastest growing 
mode and gives an order-of-magnitude estimate for the domain length by / ~ l/|<?maxl- For 
7i = J2 we obtain for the phase separation domain length 



where > A11A22 [where A;j is defined in Eq. ([8])] and the latter equality is valid for 
-^11 = A22. The domain length from Eq. (|46] ) depends on the lattice height that modifies 
the dispersion relation. In deep lattices and with stronger interactions the domain length 
becomes shorter according to Eq. (|46l i. as also numerically demonstrated in Ref. [|39l . Note 
that reaching the ground state of maximally phase-separated state would require a very slow 
transition to the unstable regime so that only the lowest energy unstable mode is seeded in the 
process. 

In Fig. [13] we show the classical field-theory estimates for the domain lengths obtained 
from Eq. (|46|) with / ~ 1 /l^maxl for different values of the nonlinearity as a function of the final 
lattice depth. The domain length rapidly decreases as the lattice becomes deeper, eventually 
saturating due to the finite range of available q values in the lowest energy band. 

In order to study the effect of quantum fluctuations of the atoms on the domain formation 
we vary the strength of quantum fluctuations in the simulations. We may continuously 
interpolate from the regime of strong quantum fluctuations to the classical mean-field limit 
by keeping the nonlinear interaction strengths Xij constant, but varying the atom number 
ll64l . This is done by changing the parameter Q.j_ajj/{Nd) for constant In the Bogoliubov 
approximation the nonlinearities Xij fix the number of excited-state atoms depleted from the 
ground state, so varying the total atom number changes the depleted fraction due to quantum 
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Figure 13. Domain length calculated from Eq. (|46] | with I ~ 1 /Ig'maxl for different values of the 
nonlinearity An - A22- The curves from top to bottom represent the values An = 2.8 x 10""', 
1.4 X 10"^, 0.7 X lO""', and 3.5 x 10""* at the lattice height s - 2. The corresponding domain 
length values for deeper lattices are obtained from the value at s = 2 by changing A,j - Xij/L 
and J in Eq. ( l46b as a function of the lattice height, as in Eqs. ( |36] | and dJTl l. The different 
curves can be considered to represent either different atom densities or different interaction 
strengths. The values of the inter-species interaction parameter y =1.5 (left) and 1.1 (right). 
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Figure 14. Quantum mechanical expectation values and uncertainties for the number of 
domains in dynamically unstable regime of the two-species condensate in phase separation. 
We vary the number of atoms per lattice site n = N/L hy changing the value of the 
parameter Q.^ajj/{Nd) while keeping the nonlinear interaction strengths Xij constant. The 
limit n — > oo corresponds to the classical mean-field result, while quantum fluctuations become 
progressively stronger as n is reduced. We find that quantum fluctuations increase the number 
of domains and the fluctuations in the number of domains. The initial value of the interspecies 
interaction parameter during the turning up of the lattice potential is in the stable regime with 
y = 0.95. After the ramping up of the potential it is changed to the unstable value of y = 1.5. 
Here L = 128 and the final lattice height i = 15. 



fluctuations. 

We show in Fig.[l4]the quantum mechanical expectation value of the number of domains 
and the corresponding quantum mechanical uncertainty, obtained by the ensemble-averaging 
stochastic phase-space simulations. The two-species mixture is miscible with y = 0.95 during 
the turning up of the lattice potential from 5, = 2 to 5 = 15, after which the system is switched 
to the unstable regime by changing the parameter value to y = 1.5. The initial nonlinearity 
Xjj = 0.60ct>R, but the number of atoms per site n = N/L in one of the species is varied from 
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Figure 15. On-site atom number fluctuations (left column) and the relative intra-species phase 
coherence between the atoms in adjacent sites (right column) for one of the atomic species in 
the unstable regime for different values of the final lattice height 5 = 10, 20, 30, 40 during and 
after the turning up of the lattice. The initial value of the inter-species interaction parameter 
during the turning up of the lattice potential is in the stable regime with y - 0.95. After the 
ramping up of the potential it is changed to the unstable value of y = 1 .5 (top row) and y -2.0 
(bottom row). 



20 to 1000, with L = 128. The limit n ^ oo corresponds to the classical mean-field limit. 
We find that both the number of domains and the fluctuations in the number of domains are 
increased due to quantum fluctuations. 

The domain boundaries between the two species may be viewed as defects and 
by introducing an electromagnetic coupling between the two components that mixes the 
atom populations provides a phase-separation scheme suitable for testing the Kibble-Zurek 
mechanism for defect formation in phase transitions [|42l l43ll . Changing the value of the 
interaction parameter through the phase transition point spontaneously breaks the symmetry 
of the system in which case the formation of the defects is expected to depend on the timescale 
of the transition, providing an interesting link to condensed matter physics and cosmology. 

3.3.2. Relative phase coherence and number fluctuations In the unstable regime there is 
a dramatic loss in the relative phase coherence between the atoms in difl^erent lattice sites 
and the growth of on-site atom number fluctuations. The enhanced on-site atom number 
fluctuations may be understood in terms of the domain formation where the spatial location 
of the domain boundaries fluctuates from one stochastic realization to another. In Fig. [15] 
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Figure 16. Stationary values of fluctuations and the density overlap between the two BEC 
components. The averaged stationary values of the relative intra-species phase coherence 
between the atoms in adjacent sites (left) and the on-site atom number fluctuations (middle) in 
the unstable regime as a function of final lattice height, obtained from the results of Fig. [15] 
The density overlap integral between the two BEC components T(f), defined by Eq. ( l26b . for 
diff'erent values of the inter-species interaction strength y. The value depends on the ratio 
between the interaction strengths ;ifii/;t'i2. 



we show the atom number fluctuations in one lattice site for one of the atomic species (left 
column) and the relative intra-species phase coherence between the atoms in adjacent sites 
(right column) for difl'erent lattice heights s = 10, 20, 30, 40 and different inter-species 
interactions y = 1.5 (top row) and y = 2.0 (bottom row). Even though there is a rapid loss 
of phase coherence due to the dynamical instability, the nonvanishing value of Ci indicates 
that even the dynamically unstable system exhibits nonvanishing phase coherence. The loss 
of phase coherence is much faster due to dynamical instability than due to the ramping of 
the lattice. The stationary averaged values as a function of the lattice height after the system 
has reached a metastable configuration are shown in Fig. [161 We also display in Fig. [16] the 
density overlap integral between the two BEC components T(t) [Eq. (l26l) 1 for different values 
of the inter- species interaction strength y that indicates the degree of spatial phase separation 
between the components. 

In Fig. [17] we show the populations in the lowest five momentum modes in the unstable 
regime. The decay of population in the ground state is fast due to the dynamical instability. 
The rapid decay of the population is a characteristic feature of the instability. We also show 
the populations of all the momentum states at a given time (t = 20000 /cor and y = 1.1) for a 
single stochastic realization and for an ensemble average over many realizations. 

4. Concluding remarks 

We demonstrated via numerical simulations how a simple approximate stochastic phase- 
space method can provide valuable information about a two-species superfluid system in 
the presence of significant quantum fluctuations. We identified contributions of quantum 
fluctuations in the pattem formation of the two-species system that results from the 
modulational instability of dynamically unstable excitations. The parameter space of the two- 
component condensate system in the lattice is especially large. Although the essential effects 
of the phase separation dynamics, the loss of the relative phase coherence between the atoms 
and the reduced atom number fluctuations were captured by the selected parameter regimes. 
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Figure 17. Occupations of the lowest five momentum states corresponding to the Fourier 
modes for different values of the final lattice height 5= 10, 20, 30, 40 in the unstable regime for 
y - 1.1, as in Fig. [15] We also show a snapshot image of the occupations of all the momentum 
components both for a single stochastic realization (middle row) and for the ensemble average 
over many realizations (bottom row). 
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the parameter space could be explored in more detail. In particular, novel phenomena would 
be observed in the case of moving condensates. The centre-of-mass motion of the atoms 
results in dissipative transport properties [l65l l64l . The stability criteria are also changed 
when the velocities of the two BECs are different and one of the BECs is not in the normal 
dispersion regime [8J. 

Moving condensates may be experimentally studied in a combined optical lattice and 
a harmonic trap by suddenly displacing the harmonic trap, e.g., by using a magnetic field 
gradient, in which case the displacement excites dipolar oscillations of atoms along the lattice 
direction [|65l . The other alternatives to create an analogous effect, e.g., are to use a moving- 
standing wave, so that the atoms experience a moving optical lattice potential ll66l introduce 
a phase shift for the hopping amplitudes of the atoms between adjacent sites, or to make the 
hopping amplitudes time-dependent by a periodically pulsating lattice. 

It would be particularly intriguing to investigate soliton structures in a two- species 
BEC. Modifying the vdXio XiilXn fT6l, e.g., by unbalanced relative populations, or inducing 
flow instabilities [fTTl [67l| would lead to the emergence of vector soli tons. The solitons 
can be persistent (metastable) or pulsating. Another potential application follows from 
the observation that repulsive interactions and quantum fluctuations lead to suppressed 
atom number fluctuations in the lattice. Two-species systems have already been used in 
experimental realizations of spin-squeezed atom interferometry [i24l and squeezing can be 
employed in interferometers to reach sub- shot-noise accuracies that are not achievable using 
classical interferometers limited by the standard quantum limit [1681169117011711 . 
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Appendix: Analytic estimates of atom number and relative phase fluctuations 

Here we provide analytic estimates for the fluctuations of the intra-species on-site atom 
number and the relative phase between the atoms in different lattice sites in the ground-state 
of an optical lattice. We introduce the atom number and phase operators for each site as 
in Ref. ||471 . In the calculations we use the classical Bogoliubov mode functions for the 
two-species BEC that are employed in the initial state decomposition of the stochastic field 
in the TWA simulations. By means of replacing the stochastic mode amplitudes a"^+ in 
Eq. (|29l ) by the annihilation and creation operators (or^^i, [Q;^"|i]^), so that 0"^^+ a^^i and 
[0"^!+]* ~^ \.ctq,±V ^ we may obtain analytic estimates for the ground-state properties, provided 
that the inter-species correlations can be ignored. 

We may derive the atom number operator at the site / of the species j by expanding 
the atom population n\^^ of the species j in the site / to first order in fluctuation terms. We 
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n^^ - ^jl^- denotes the ground state atom number per site of the species j and = d^l^-v, 
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(A.l) 



^^(6^^' + [<5>;^V) = E + H.C.) , (A.2) 



(i) 



,0-)_„U) 



We may introduce the corresponding phase operator at the site / as 



r'-;!'„a'-;:'„e"'' - H.c.j , 



(A.3) 



for which the commutator [hy, (py'] = i and we have defined r^^ = Ug'r, + Vq^lj. Then the on-site 



,0') , „o-) 



atom number fluctuations in the Ith site (An9^)^ and the relative phase fluctuations between 



the atoms in the kth and Ith site (A(;c[f respectively, read 



L ^ 2 



eq 1 2nq- + 1 2n^_+ + 1 



(A.4) 



q{k - 1) 



^ Z 27 I^^"'^'- + ^ ^^""^'^ + ^^"^1 '''' 



(2n,,, + 1) 
q{k - I) 



(A.5) 



where + is the thermal population of the phonon mode {q, +) in the lattice given in Eq. (I32l) . 
In the second line of Eqs. (IA.4I) and (IA.5I) we have used the specific relations for m^-[],, v^], 
and given by Eqs. (I22l) . (I23l) and (|2T]) . respectively. These were obtained by choosing a 
particular set of parameter values {^u = X22, Ji = J2 = J, M = N2, etc.). The single-particle 
energy = 4ysin^(^/2). 

At r = we have n^^r] = and we can evaluate the expressions analytically by replacing 
in the continuum limit the momentum sums by integrals 



For the atom number fluctuations we obtain 

n^j^ 



(Any'y = — [arctan (i+) + arctan ( A.)] 



n 



(A.6) 



where (An = A22 > A^ > 0; A;y is defined in Eq. ®) 
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Figure Al. The analytic estimates of the ground-state fluctuations. On left, the scaled on- 
site atom number fluctuations n'q'j, = (A«J-'V yjni in one of the lattice sites, obtained from 
Eq. ( |A.6t . In the middle, the relative phase coherence between the atoms in the adjacent sites 
C'j, calculated from Eq. (lA.llt . In the both cases the upper curve corresponds to y = 0.9 and 
the lower curve to 7 = 0.1. On right, we show the phase coherence C'j along the lattice as a 
function of the lattice site separation j, calculated from Eq. (lA.lll i. The curves from the top 
represent (s = 20, 7 = 0.1), (s = 30, y = 0.1), (s = 40, 7 = 0.9) and (s = 40, 7 = 0.1). 



In the nonlinear limit of (An + A12) » J this simplifies to 

(A„W)2 ^ ^ + (A.8) 

n \ VAii + A12 VAii - A12/ 
Similarly, for the relative phase fluctuations between the atoms in the adjacent sites {k-l = 1) 
we obtain 

1 



hrA^.^') "''"^'^^ ^[h') "'^'"^'-^ 



.(A.9) 



In the nonlinear case (An ± A12) » J this reads 

= ( .) ^ ( VA11+A12 + VA11-A12) . (A.IO) 

Similar relationships as Eq. (|A.9I ) may also be calculated for difl^erent values of k in (A^|^^V 
in the continuum limit. 

The analytic expression (IA.91 ) for the relative phase fluctuations can be used to evaluate 

C; = <exp [z (^^+, - (p^\) ^ exp [- ((^^+, - (p^f^ /2] , (A.l 1) 

as displayed in Fig. lAll in which case we show the relative phase coherence between the atoms 
in the adjacent sites C\ together with the on-site atom number fluctuations. We also use the 
continuum limit approximation to calculate {iSpk^j - 'Pkf'^ for different values of j in order to 
obtain the coherence along the lattice C'j as a function of the site separation j. The parameters 
of Fig. lAll are the same as those used in in the TWA simulations for the dynamically stable 
r = cases, with the dependence of the hopping amplitude J and the nonlinearity n = xii 
on the lattice height determined by Eqs. (|36l ) and (|37] ). respectively, where we set the value 
Xn = O.^ojr at 5 = 2. The atom number A'^ = 2560 and the number of sites L = 64. 

In Fig. lAll the relative phase coherence between the atoms in the adjacent sites decreases 
rapidly as a function of the lattice depth and increases as the inter-species interaction strength 
7 is increased closer to the onset of the phase-separation instability. The numerical values 
of the nearest-neighbour coherence are very close to those of the TWA simulations in Fig.[Tl 
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but the long-range coherence values are higher than in the TWA case [Fig. [21 . One should 
note, however, that C[ does not include the atom number contributions incorporated in the 
definition of C\ [Eq. (|4T)) ] which is used in analyzing the relative phase coherence in the TWA 
numerics. 

We can implement a nonlinear least square fit for the coherence along the lattice C'j using 
a trial function 

= ai exp(-a2x'") + a^ , (A. 12) 

and determine the coefficients a,. An accurate fit for the calculated values in Fig. lAll 
with y = 0.1 is obtained for (ai,a2,ai„a4) = (0.11,0.55,0.16,0.92) for s = 20 and 
(0.24,0.44,0.38,0.77) for s = 40. These indicate asymptotic values for the coherence for 
large spatial separations 0.92 and 0.77, respectively. 

The on-site atom number fluctuations from the analytic estimates in Fig. lAll are also very 
close to the numerical TWA simulation results of Figs. |5]and |6l The on-site atom number 
fluctuations are enhanced as the inter-species interaction strength increases. 
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